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Abstract 

In an optical trap, the ground state of spin-1 Bosons such as 23 Na, 39 K, and 
87 Rb can be either a ferromagnetic or a "polar" state, depending on the scat- 
tering lengths in different angular momentum channel. The collective modes 
of these states have very different spin character and spatial distributions. 
While ordinary vortices are stable in the polar state, only those with unit 
circulation are stable in the ferromagnetic state. The ferromagnetic state also 
has coreless (or Skyrmion) vortices like those of superfluid 3 He-A. 
Two months ago, the MIT group has succeeded in trapping a 23 Na Bose condensate by 
purely optical means [|]]. This experiment has opened up a new direction in the study of 
confined dilute atomic gases. In conventional magnetic traps, the spins of the alkali atoms 
are frozen. As a result, even though the alkali atoms carry spins, they behave like scalar 
particles. In contrast, the spin of the alkali atoms are essentially free in an optical trap. The 
spinor nature of alkali Bose condensate can therefore be manifested 0. Because of the wide 
range of hyperfine spins of the alkali Bosons and Fermions, the optical trap has provided 
great opportunities to study dilute quantum gases of atoms with large spins. 

The purpose of this paper is to point out the general properties of the spinor Bose 
condensates. As we shall see, they possess a whole host of quantum phenomena absent 
in scalar condensates. These include changes in ground state structures with interaction 
parameters, vector and quadrupolar spin wave modes, topological and energetic instability 
of doubly quantized singular vortices, and the existence of coreless (or Skyrmion) vortices. 
All these results are simple consequences of the effective low energy Hamiltonian of the 



system. The derivation of this Hamiltonian (which applies to Bosons and Fermions with 
arbitrary spins) and the realization of its limitations are therefore crucial for our discussions. 

For simplicity, we shall consider Bosons with hyperfme spin / = 1. This includes alkalis 
with nuclear spin I = 3/2 such as 23 Na, 39 K, and 87 Rb. Alkali Bosons with / > 1 such 
as 85 Rb (with / = 5/2) and 133 Cs (with / = 7/2) have even richer structures and will be 
discussed elsewhere. To illustrate the fully degenerate spinor nature, and as a first step, 
we shall only consider the case of zero magnetic field. To reach a good approximation of 
this limit, the Zeeman energy must be much smaller than the interaction energy at the 
center of the trap (which is essentially the chemical potential fi) . The condition is therefore 
7 5 << /i, where 7 is the gyromagnetic ratio, and B is the magnetic field. For clouds with 
10 6 atoms, with scattering length a sc = IOOA and trap frequency 200Hz, this condition gives 
B « 10~ 3 gauss. This bound can be further increased by increasing the trap frequency. 
Since the current capability of magnetic shielding can reach 10~ 5 gauss, a good approximate 
of the zero field limit is attainable. 

Our key results are : (1) The spinor ground state depends crucially on the s-wave 
scattering lengths {op} of various elastic scattering channel with total hyperfme spin F. 
When a > a 2 , the ground state is a ferromagnet denoted by a normalized spinor ( + = 
(1, 0, 0). When a 2 > a , the ground state is a non- magnetic state "polar " state denoted by 
a normalized spinor ( + = (0, 1,0). Current estimates of a and a 2 indicate that the spinor 
condensates of 23 Na and 87 Rb may be a polar and a ferromagnetic state respectively. 
(2) The polar state has a density mode (ujq) and two degenerate "vector- like" spin wave 
modes (w±), all of which are of Stringari form ||. For large clouds, the density mode is 
universal (i.e. independent of interaction), whereas uo± = ^-^^o- The ferromagnetic state 
has a density mode similar to that of the polar state. It also has a "vector-like" and a 
"quadrupolar-like" spin wave mode, corresponding spin flips of one and two units. Both 
spin waves modes obey a Schrodinger equation and reside near the surface of the cloud. 
The quadrupolar modes is separated from the vector mode by a gap of the order of the trap 
frequency. 
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(3) In case of ferromagnetic states, vortices with even circulation (±2mn, m > 1) are 
neither topological or energetically stable. Only vortices with unit circulation (±27r) have 
topological stability, a phenomena identical to that of superffuid 3 He-A. 

Optical Trap Potential and Their Spin Dependence: We begin by examining the confining 
potential produced by a laser with a linearly polarized electric field £(£)x and frequency u. 
The potential seen by an atom in hyperfine state |l,m > is U m = — |a m £ 2 , 

where (..) denotes the time average, a m is the polarizability of the state \l,m >, hug m > 
is the excitation energy from the ground state |l,m > to excited states \£ >, eXe m is the 
dipole matrix element between |1, m > and \t >, and Z is the atomic number. We have also 
written down the dipole sum rule in eq.fll]). For a red detuned laser, to < ue m , the atom will 
be attracted to region of strong field. A focused laser beam will therefore form a harmonic 
trapping potential near the point of maximum intensity. 

In general, the spin dependence of the trapping potential is weak. The reason is that 
to one part in 10 2 to 10 3 , the dipole sum rule is saturated by the nS to nP transitions [|J, 
where nS is the ground state electronic configuration. The sum rule in eq. (|l|) can therefore 
be approximated as Y^eep\-^im\ 2uj tm — of?-^- I n addition, the excitation energies uj£ m for 
different nS to nP transitions differ from each other only by the fine splitting, which is 
typically 10~ 3 smaller than the energy difference (hups) between nS to nP levels before 
they fine split. Thus, for a far detuned laser, eq.(|l|) can be approximated as 




where Afi ne is the fine splitting and Vt = oops — lo is the detuning frequency. To the extent 
that the transitions nS — > n'P (n' > n) and the ratio Af ine /Q can be ignored, the sum in 
the first term in Eq.(|2|) becomes a constant and the correction in the second term can be 
set to zero. This implies that a m , hence the trapping potential U m , is independent of m. 

Effective low energy Hamiltonian : Alkali atoms (Bosons or Fermions) have two hyperfine 
multiplets (denoted as jhigh and fi ow ), with energy difference many orders of magnitude 
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larger than the trap frequency. Since the interaction between two alkali atoms depends on 
their electron spins (singlet or triplet), the hyperfme states of the atoms can be changed after 
the scattering. However, when the system is very cold, two atoms in fi ow will remain in the 
same multiplet after the scattering since there is not enough energy to promote either atom 
to fhigh- In contrast, energy conservation does not limit the production of fi ow states from 
the scattering of two fhigh states. As a result, in an optical trap, all atoms in the ground 
state will be in the lower multiplet. The low energy dynamics of the system will then be 
described by a pairwise interaction (between atoms 1 and 2) that preserves the hyperfme 
spin of the individual atoms and is rotationally invariant in the hyperfme spin space. The 
general form of this interaction is 

2/ 

V{v 1 -v 2 ) = 5{v 1 -v 2 )Y j 9fV F} g F = Anh 2 a F /M (3) 

where M is the mass of the atom, V F is the projection operator which projects the pair 1 
and 2 into a total hyperfme spin F state, and a F the s-wave scattering length in the total 
spin F channel. For Bosons (or Fermions), symmetry implies that only even (or odd) F 
terms appear in eq.(|3]). 

For a system of / = 1 Bosons, we have V = giV 2 + goPo- Likewise, the relation 
Fi ■ F 2 =Ef / =o Wf, Xf = \ [F(F + 1) - 2/ (/ + 1)] becomes F x • F 2 = V 2 - 2V . The 
relation J^f'PfLo = 1 becomes 1 = V 2 + Vq. We then have (dropping the 5-function in 
eq.(D), 

V = c + c 2 ¥ l ■ F 2 , (4) 

where c = (go + 2g 2 )/3 and c 2 = (g 2 — go)/ 3. The Hamiltonian in the second quantized 
form is then 

H = I dr Wa+ ' Wa + + + f^VtfFa* • F a , 6 ,#V&) (5) 

where ip a { r ) is the field annihilation operator for an atom in hyperfme state \l,a > at 
point r, (a = 1,0, —1), and U is the trapping potential. It is straightforward to generalize 
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eq.(||) to higher spins. By noting that (Si • 82)™ = Y?J=o ^f'Pf, eq.(|J) can be written as 
V = YJn=o c n (Fi ■ F 2 )™, where c n is a linear combination of the gp's. Alternatively, it is useful 
to express the projection operator Vf in the second quantized form as Vf = J2a=-F OpaOpa, 
where C>Fa= J2 ai ,a 2 < Fa\f,ai,f,a 2 > ip ai ip a2 , and < Fm\f ,ai; f, 0,2 > is the Clebsch- 
Gordon coefficient for forming a total spin F state from two spin-/ particles. The pair 
potential is then V = J2f=o9f J2a=-F OpjOpa- Only even and odd F terms appear in V in 
case of Bosons and Fermions. 

Ground State Structure : It is convenient to write the Bose condensate ^f a { T ) =< ^a( r ) > 
as # a (r) = 



'n(r)( a (r) } where n(r) is the density, and ( a is a normalized spinor ( + ■ ( = 1. 
The ground state structure of ^a(r) is determined by minimizing the energy with fixed 
particle number, i.e. 8K = 0, K = 5 < H — fj,N >, where fi is the chemical potential, 



K 



( ti 2 h 2 n 2 r 

dr 2M (Vv ^ )2 + 2M (VC)2n ~ [/i ~ U{r)] n + Yl Co + C2<F> ' 



(6) 



and < F >= £*F &(&. It is obvious that all spinors related to each other by gauge transfor- 
mation e ld and spin rotations U(a, (3, r)=e~ lFz0t e~ tFvt3 e~ lFzT are degenerate, where (a,f3,r) 
are the Euler angles. There are two distinct cases : 

I : Polar state : This state emerges when c 2 > (i.e. g 2 > go). The energy is minimized by 
< F >= 0. The spinor ( and the density n° in the ground state are given by 



C = e m U 



je 



^ — -j^e iQ sin/? ^ 
cos/3 



n°(r) = —Ui-U(t)-W(t)) 

Co 



(7) 



where W(r) = ^ V ^T - Note that ( is independent of the Euler angle r. The symmetry 
group of polar state is therefore U{\) x S 2 , where U(l) denotes the phase angle 9, and S 2 
is a surface of a unit sphere denoting all orientations (a,/3) of the spin quantization axis. 
II : Ferromagnetic state : This state emerges when c 2 < 0, or go > g 2 . The energy is 
minimized by making < F > 2 = 1. The ground state spinor and density are 
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C = e a U 
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V 



e ia sin 2 § 



n°(r) = — (u- U(r) - W(r)) 
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The direction of the spin is < F >= cosaz+ sina(cos/5x + sin/3y). The combination (6 — 
7) in eq.(P) clearly displays a "spin-gauge" symmetry [0], i.e. the equivalence between 
phase change and spin rotation. Because of this symmetry, the distinct configurations of ( 
(including the gauge) are given by the full range of the Euler angles. The symmetry group 
is therefore SO (3). As we shall see, this difference in symmetry between the polar and the 
ferromagnetic state leads to a fundamental difference in their vortices. 

According to the latest estimate of J. Burke, J. Bohn, and C. Greene ||, the scattering 
of 23 Na are 02 = (52±5)a# and ao = (46±5)a#; and those for for 87 Rb are 02 = (107±4)a# 
and ao = (110 ± 4)a#, where a# is the Bohr radius. Because of the overlapping error bars 
of (Z2 and ao in each case, one can not say for sure about the nature of their ground states. 
However, if the inequalities suggested by current estimate (a 2 > ao for 23 Na and a > a 2 for 
87 Rb) are true, then the condensates of 23 Na and 87 Rb are polar state and ferromagnetic 
state respectively. 

Collective modes of trapped spinor Bose condensates: The equation of motion in zero 
field is 



i%d t ip r , 



h 2 



2M 



V 2 ^ m + [U(r) - n\$ m + c (^+4) $m + c 2 (^+F o6 4) • (F^) f 



(9) 



To study the elementary excitations, we write i[j m = ^>° m + 4> m and linearizing eq.(|9]) about 
the ground state 

I : Polar State : Without loss of generality, we take z as the spin quantization axis and 
( T = (0, 1,0), with the subscript "T" denotes the transpose. Using the expression of n° in 
eq.(|7|) and the fact that < F >= 0, eq.([5]) becomes 
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To linear order in m , the density and spin fluctuations (Sh and 6M±) are related to 0o and 

andM_ = M±. 



0± as 5n(r) = ^/n°(r) (0 O + 5M+ ee S(M x +iM y ) = ^/n°(r) (fa + 
Denoting the frequencies of fa and 0-t as u and u;±, it is easy to see that they are all of 
the Bogobuibov form in the homogenous case (W = 0), where Kloq = we k (ek + 2con°), 



u± = Je k (e k + 2c 2 n°), and e k = h 2 k 2 /(2M). 

In a harmonic trap, eq.(|10"D is identical to the equation of the collective mode of a 
scalar Boson. Using the method of Stringari ||, it is straightforward to show that for large 
cloud (where n° in eq.(^) is well approximated by the Thomas Fermi expression n°(r) = 
[ji — U{y)]/cq by ignoring W(r) 0), eqs. fllOl) and (PI) can be written as 

d 2 5h = V (c n°V5n) , <9 2 <5M± = V (c 2 n°V5M±) . (12) 

Stringari has shown that the density mode in eq. (|T2"D has a universal spectrum (i.e. in- 
teraction independent) with power law (hence extended) wavefunctions ||. Since the spin 
waves modes 5M± obey exactly the same equation as 5n except that c is replaced by c 2 , 
the quantum numbers and the wave functions of the spin wave modes are identical to those 
of the density mode, and 



c 2 a 2 — a o 

u± = —uj = - ■ to . 

c 2a 2 + a 



(13) 



II : Ferromagnetic State : We shall take ( T = (1,0,0). Using eq.(^) and the fact that 
< F >= z, eq.(P) becomes 
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To linear order in a , the density, spin, and "quadrupolar" spin fluctuations are given by 
5h = y^n°(fa + (frf), 5M_ = y/npipQ , and 5M 2 = 2v / n°0Ti- The frequencies of these modes 



will be denoted as u>i, uo, and U-i respectively. In the homogenous case (W = 0), the density 
mode has a Bogoluibov spectrum hu = ^J^+2g^n°e^, ek = {Tik) 2 /(2M). The spin wave 
SM_ has a free particle spectrum uj = e^. The frequency of is free particle like with 
a gap, uj Q = e k + 2\c 2 \n°. 

In a harmonic trap, the density mode 5n assumes the Stringari form because 0i obeys the 
same equation as scalar Bosons. The spectra of both O and 0_i are given by Schrodinger 
equations with potentials ^(r) and W(r) + 2|c2|n°(r) respectively. Outside the cloud, the 
potentials for both 0o and 0_i reduce to the harmonic potential U(r) — \i. For a large cloud, 
n° is well approximated by the Thomas Fermi expression n° = \ii — XJ\j g 2 except near the 
surface of the cloud. Calculating W with n°, one finds that it becomes more attractive 
near the surface. The low energy modes of 5M_ are therefore confined near the surface. 
For the quadrupolar spin waves, 5M 2 , the low energy modes are even more confined to the 
surface because of the additional potential 2|c2|n°(r), which is much more repulsive than 
W inside the cloud. The spectrum of the surface modes of 5M_ and therefore mimic 
their homogeneous counterpart, with the energy of SM^_ shifted up from that of 5M- by an 
amount ~ \c2\n ~ (i^l/fl^)/^ = (\c2\/92)(R/o,t) 2 ^t, where R is the size of the cloud, and 
ax = \Jh~/ Mojt is the trap length. 

Intrinsic stability of singular vortices with circulation (£ > 1) in the ferromagnetic state: 
The fundamental difference between the vortices of the polar and the ferromagnetic states 
can be illustrated by their superffuid velocities, v s = From eqs.(0) and (H), they 

are 

(v s ) polar = (v,) /erro = ^(V(^-r)-cos/5Va), (15) 

where we have assumed that the Euler angles (a, (3, r) and 9 are spatially varying functions. 
Unlike the polar state, the superffuid velocity v s of the ferromagnetic state depends on spin 
rotations. This leads to the following remarkable property of the Bose ferromagnets : if 
too much vortex energy is stored in one spin component, the system can get rid of it by 
spin rotation. To illustrate this phenomenon, consider the following family of spinor states 

8 



{ty a (t) = Vn°( a {t)} parametrized by a parameter t between and 1, 

e«) = ( e «"W (f ) ,e-*V2si„ (f ) cos (f ) ,W (f )) , (16) 

where m > is an integer, is the azimuthal angle, and n° is the equilibrium density for 
the vortex state ( T (t = 0) = (e 2m<?il , 0, 0) with £ = 2m circulation. As t evolves from to 1, 
this 2m7r-vortex evolves continuously to the vortex free state ( T (t = 1) = (0, 0, 1) with a 
spin texture < F >=cos(-7rt)z + sin(7rt)[cos(m0)x + sin(m0)y]. 

Eq.(|T6|) shows that the 2m<f) vortex is topologically unstable. It can only be stabilized 
by energy barriers preventing its collapse. Such barriers, however, are non-existent. From 
eq.(|5]), it is easily shown that (with c = cos (y) , s 2 = 1 — c 2 ), 

dK it h 2 



dt 2 2M 



J dm°(r)|V2m0| 2 |(c+ 2c 3 )s < 0, (17) 



for < t < 1. The fact that K decreases monotonically with increasing t shows the absence 
of energy barriers. 

By multiplying eq.(^) by e*^, one also obtains a family connecting a vortex with 2m + 1 
circulation ( T (t = 0) = e^ 2m+1 ^(l, 0, 0) to a vortex with unit circulation ( T (t = 1) = 
e l( ^(l, 0, 0). It can be easily shown as in eq.(|TTD that dK/dt < for m > 0. This means that 
all (2m + l)7r-vortices will collapse into a 27r-vortex. The 27r-vortex, however, can not be 
deformed continuously into a uniform state ||. 

Coreless (or Skyrmion) vortices: Eq . (|I5|) shows that spin variations in the ferromag- 
netic states in general leads to superfiows. To illustrate further, consider the condensate 
C(r) T =(cos 2 |, y/2e l sin ^cos^, e 2l ^sin 2 ^) where (3 = (3(r) is an increasing function of r start- 
ing from p = at r = 0. The spin texture and superfluid velocity of this condensate are 
both cylindrically symmetric, < F >=zcos/5 +sin/9(cos0x-|-sin0y), and v s = — cos/3) 0. 
When /3(r) reaches tt/2, v s becomes the velocity field of a singular vortex. However, the 
singularity of the usual 2tt vortex is absent because v s vanishes instead of diverges at r = 0. 
This phenomena is identical to that of superfluid 3 He-A, which has exactly the same angular 
momentum texture, superfluid velocity, and topological instability 11 ||10||. In the case of 



3 He-A, it is known that external rotations can distort the texture so as to generate a velocity 
field v s to mimic the external rotation as closely as possible 0. Such textural distortions 
will occur here for the same energetic reason, physical reason. 

We have thus established Statements (1) to (3). In current experiments, the condensates 
are first produced in a magnetic trap and then loaded into an optical trap. Because of the 
way they are produced, the condensates will in general carry a net magnetization. If the 
spin relaxation time is sufficiently long, the condensate will behave like a ferromagnetic state 
even its ground state is a polar state. In this case, the vortex phenomena discussed above 
will apply. 

(Notes added: At the time of submission of this paper, a preprint by T. Ohmi and 
K. Machida has just appeared (|cond-mat / 9803 1601) ■ These authors have studied the same 
model in finite field and in the absence of a trap. Our zero field results of homogeneous 
systems agree with each other. ) 
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